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ABSTRACT 

Context. We have conducted a 3D MHD simulation of the solar corona above an active region in full scale and high resolution, which 
shows coronal loops, and plasma flows within them, similar to observations. 

Aims. We want to find the connection between the photospheric energy input by field-line braiding with the coronal energy conversion 
by Ohmic dissipation of induced currents. 

Methods. To this end we compare the coronal energy input and dissipation within our simulation domain above different fields of 
view, e.g. for a small loops system in the active region (AR) core. We also choose an ensemble of field lines to compare, e.g., the 
magnetic energy input to the heating per particle along these field lines. 

Results. We find an enhanced Ohmic dissipation of currents in the corona above areas that also have enhanced upwards-directed 
Poynting flux. These regions coincide with the regions where hot coronal loops within the AR core are observed. The coronal density 
plays a role in estimating the coronal temperature due to the generated heat input. A minimum flux density of about 200 Gauss is 
needed in the photosphere to heat a field line to coronal temperatures of about 1 MK. 

Conclusions. This suggests that the field-line braiding mechanism provides the coronal energy input and that the Ohmic dissipation 
of induced currents dominates the coronal heating mechanism. 

Key words. Sun: corona — Magnetohydrodynamics (MHD) — Methods: numerical — Sun: UV radiation 


1. Introduction 


The coronal heating mechanism is poorly understood, although 
various theories exist - many of which are listed in the review of 
Klimchuk] ( |2006| ). A coronal heating mechanism including the 
actual dissipation is field-line braiding due to horizontal shuf¬ 
fling motions in the photosphere that lead to induced currents in 
the corona where these currents heat the plasma though Ohmic 
dissipation ( [Parker 1972| l. Models of solar active regions indi¬ 
cate that this mechanism can produce a loop-dominated corona 
( jGudiksen & Nordlund|2002l|2005b|aj l. 

Statistical analyses of synthesized spectra from these mod¬ 
els match statistical properties of observed spectra ( [Peter et ak] 
2004 [2006 1. The Doppler shifts in these models were able to 
reproduce the observed persistent redshifts in the transition re¬ 
gion, but failed to show the observed blueshifts at higher temper¬ 
atures above 500000 K ( Peter|1999l Peter & Judge|1999]l. This 
is also true for more recent active region models ( [Zacharias et al.[ 
|2011b|a[ ). Though quiet Sun (QS) models do show blueshifts at 
high temperatures, the magnitudes of their red- and blueshifts 
are too small ( [Hansteen etal.|20l0) l. We present here a model of 
an active region in a full-scale computational domain that imple¬ 
ments the field-line braiding process and that is able to repro¬ 
duce the redshifts in the transition region and the blueshifts in 
the corona. 


We can follow the magnetic field through the model corona 
and extract any quantity along that field line, which is otherwise 
inaccessible by coronal observations. We find that a minimum 
magnetic flux density of 200 Gauss is necessary at the footpoints 


of these loops to produce a sufficient upwards directed Poynting 
flux to heat the coronal part of these loops by Ohmic dissipation. 


2. Coronal model and properties of energy input 

In our study ( [Bourdin et al.[201^ , we have conducted a 3D mag- 
netohydrodynamical (MHD) numerical simulation of the solar 
corona above an active region (AR). The model includes grav- 


ity, radiative losses following 

Cook etal. (|l989]l, tield-aligned 

heat conduction (Spitzerj 

1962 

1 , and a uniform magnetic diffu- 

sivity rj = lO'*^ m^/s. T 

he plasma pressure and the magnetic 


fleld in the corona are described self-consistently. In order to do 
this, it is important to set up a realistic energy balance, including 
necessary coronal energy sinks. The Spitzer-type heat conduc¬ 
tion transports energy downwards to the chromosphere and also 
powers evaporation of colder plasma there. 

We drive the model from below by photospheric obser¬ 
vations, using Hinode/SOT-NFI ( [Kosugi etal.[[2007| [Tsuneta| 
et al. [2008 1 line-of-sight magnetograms and horizontal veloci¬ 
ties. Large-scale motions of magnetic patches and granulation¬ 
like small-scale velocities introduce magnetic stress at the mag¬ 
netic held footpoints. These perturbations propagate along the 
held from the photosphere into the corona and carry magnetic 
energy as Poynting flux. Ohmic dissipation of induced electric 
currents turned out to be the main contribution to the model 
corona heating mechanism. 

With the atomic database Chianti (Dere et al. [1997[ Young 


etal. 


etal. 


2003 [) we deduce the synthetic emission following Peter 


( 2004[ [2006[ l and compare it to real observations of the 
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Fig. 1. Active region observed by the Hinode satellite on 14 Nov 2007. The left panel shows a line-of-sight magnetogram (saturation level: ±300 G) 
of the AR including magnetic field lines along the most prominent coronal loops (CL 1-2 and SL 1-3, see |Bourdin etal.|2013^ . The middle panel 
shows the Fe xiii EUV emission synthesized from the model data. The footpoints of the loops are marked with circles, where blue indicates plasma 
upfiows and red downflows. In right panel we display the vertical component of the Poynting flux at 3 Mm height. Upwards oriented vertical flux 
is colored in blue and downwards flux in red, both in linear scale with a saturation level of ±50000 W/m^. The solid square indicates the hot AR 
core area, while the dotted rectangle outlines the FOV of only the short loops SL 1-3, see also Sect.|3.4| 


same active region, like the EUV and X-ray emission spectra 
from highly ionized atoms and their corresponding Doppler line- 
shifts, taken by Hinode/EIS ( Culhane etal.|[2007] l. We found a 
substantial match of coronal stractures and their flow dynamics 
between the simulated and observed loops system. 


2.1. Boundary conditions 

Our model is periodic in the horizontal directions. The lower and 
upper boundaries of the computational domain are closed to any 
plasma inflow or outflow, and the temperature gradient is zero at 
the top and bottom boundaries to avoid any artificial energy flow 
into or out of the domain. The density is set according to the 
temperature, so that the boundary of the physical domain is in 
hydrostatic equilibrium. We extrapolate the magnetic field at the 
upper boundary with a potential field into the three ghost layers. 

For the magnetic field at the lower boundary, we use the 
potential-field extrapolation into the interior of the Sun, which is 
in principle identical to a potential-field extrapolation to the ex¬ 
terior, except that we do not decrease the contrast exponentially 
depending on the distance from the physical boundary layer; in¬ 
stead, we increase the contrast of the vertical magnetic field with 
depth. The pressure scale height below the photosphere is about 
300 km, but for numerical reasons this scale height is stretched in 
our boundary condition to about 1000 km. On the other hand, for 
the relatively short distance of about 300 km that we need to ex¬ 
trapolate into the solar interior, this choice seems to be relatively 
safe because the typical length-scales of advectional structures 
in the upper convection zone are usually larger than 1000 km in 
the vertical direction. Furthermore, a small error in the horizon¬ 
tal field component at the photospheric level would have almost 
no influence on our model corona. 


2.2. Initiai condition and switching on 


We start with a stratified atmosphere that we combined from 
temperature and density profiles of the solar interior (|Stix 


2004|l, the chromosphere (|Fontenla etal.||1993]l, and the corona 


(November & Koutchmy| 19961. We relax this combined profile 


into equilibrium using a ID hydrodynamic model of a vertical 
atmospheric column together with a velocity damping and the 
same diffusion parameters as we use for the 3D MHD model. 

Granulation has a typical life time of about 5 minutes and the 
Alfven crossing time of the magnetic energy input from the pho¬ 
tospheric driving to reach the corona is about 30 minutes. From 
minute 10 to 35, we smoothly switch on the radiative losses and 
the Spitzer-type heat conduction, to compensate for a lack of 
coronal heating during this initial phase. All physical processes 
of the model are fully active after minute 35. To damp switch-on 
effects in the 3D MHD model, we use a velocity damping that 
fades out smoothly until minute 40. 

After the Alfven crossing time plus about three granular life 
times, we see a strong rise in the coronal Ohmic heating around 
minute 48 in the simulation. At that time, the maximum temper¬ 
ature from the initial atmosphere has dropped from 1.2 MK to 
about 0.4 MK and then rises again to about 1.4 MK until minute 
60, where it reaches a kind of plateau. The small loops system 
(SL 1-3, see Sect. O' is developed at minute 63 and was then 
stable in the model for about 7 solar minutes. 


3. Energy input viewed at iarge 

3.1. Overview of whoie box and subvoiumes 

The AR we investigate in this work has a size of about 200 x 
100Mm^. Observations exist from the photosphere, in particu¬ 
lar LOS magnetograms and calibrated vector magnetograms of 
a smaller FOV; coronal observations were also made co-spatial 
and co-temporal in the range of X-rays and extreme ultra violet 
(EUV), see Fig.[2 For the simulation we use the full FOV shown 
in Fig.[T] including some surrounding QS, so that effects of the 
horizontal periodicity within our model are minimized. Mainly, 
the FOV can be divided into three different regions: the full FOV 
including the QS, the AR core area, and the small loops system 
within the AR core, see Fig.[2 While restricting our analyses to 
the small loops system (SL 1-3) for example, we can study a re¬ 
gion with enhanced coronal energy dissipation and strong EUV 
emission. 
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Fig. 2. Horizontal averages of the volumetric Ohmic (red) and viscous 
heating (blue) plotted as a height profile, where the line styles (solid 
and dashed) depict the different FOVs given in brackets in the leg¬ 
end. The dash-dotted lines basically describe the same physical quantity 
Wtoremz = u • (V X B X B) /yu, while the Lorentz work (green) and the in¬ 
duction work (orange) are of opposite sign. Due to energy conservation, 
the sum of all these quantities must be equal to the Poynting flux diver¬ 
gence, e.g. for the SL 1-3 FOV (black dotted). 


Fig. 3. Similar to Fig.[^ but for the heating per particle. Horizontal 
averages of Ohmic (red) and viscous heating per particle (blue) versus 
height. The dotted curves indicate the maximum values at each height, 
while the solid lines show the corresponding mean value. Overplotted is 
the height profile of the maximum temperature (green dash-dotted). The 
different FOVs are symbolized as line styles and are given in brackets 
in the legend. 


3.2. Contribution to the internet energy 

To investigate the spatial distribution of the energy input into 
the model corona, we look at the individual heating terms in a 
horizontal average at a given height. The driving in the bottom 
boundary layer provides the energy input by Poynting flux into 
the model. 

In Fig. 1^ we show which energy dissipation mechanism is 
active at which height for different fields of view (FOV). While 
below the corona all profiles lie relatively close together, looking 
at the full FOV data (dashed lines) we see that in the corona the 
Ohmic heating dominates over the viscous heating at almost all 
heights. In general, the coronal part of the profiles drops roughly 
exponentially with height. This applies for the box average as 
well as for individual loops (SL 1-3). The three loops SL 1-3 
show a very similar variation of the heat input along the held 
lines (and thus with height), which is why we show here only one 
curve averaging the three loops. The result of the energy input 
dropping with height for individual held lines is consistent with 
I van Wettum etal.| ( |2013| l who found that the volumetric heating 
drops exponentially with height also for individual held lines. 

The scale height of the Ohmic and viscous heating (aver¬ 
aged over the full box FOV) are very similar and uniform in the 
whole corona. When we restrict the FOV to the hot AR core 
loops SL 1-3 (solid lines) the Ohmic heating profile gets more 
structured. Between the heights of 5 and 30 Mm, where the hot 
loops are located, we see a significant increase of the Ohmic 
heating compared to the full FOV. 

In addition, we see a shortage of the Ohmic heating above 
the hot AR core loops in heights between 40 Mm and 90 Mm in 
relation to the full FOV average (cf. Fig.[^. This can be under¬ 
stood by assuming that the energy transport is directed mainly 
upwards and that an increased dissipation in a lower layer would 
decrease the possible energy input to an upper layer. Of course 
there is more to consider; for example, most of the held lines that 
close well above the hot AR core subvolume are rooted outside 
that area, so that energy can be transported there while following 
the magnetic held. 


Both the full and reduced FOV profiles in Fig.|^ approxi¬ 
mately join up again at heights above 90 Mm, which indicates 
that either some energy is transported to above the hot AR core 
loops subvolume from outside. Alternatively, this upper part of 
the corona could simply not yet have become fully independent 
of the initial condition because the Alfven travel time for mag¬ 
netic perturbations from the bottom layer into the upper corona 
is significantly longer than into the hot AR core loops system. 

The profile of the divergence of the Poynting flux mainly 
follows the sum of the Ohmic heating and the Lorentz work (see 
Fig.|^. Between 10 and 20 Mm we And a region that is in aver¬ 
age dominated by induction work. For the full FOV average, this 
induction region would expand farther from 10 to about 55 Mm. 

The viscous heating also shows deviations between the two 
FOV averages, but the difference is much smaller than the devi¬ 
ations seen in the Ohmic heating (red and blue lines in Fig.|^. 
Also for the small FOV, the viscous heating is weaker between 
heights of about 10 to 100 Mm than for the full FOV. Above 
6 Mm height the viscous heating in the small FOV (containing 
SL 1-3) is always at least a factor of 5 weaker than the Ohmic 
heating or the Lorentz work. Therefore, the viscous heating only 
provides a minor contribution to dissipating the Poynting flux in 
the corona. Altogether, this shows that the heat input to the hot 
AR core loops in our model is largely dominated by the Ohmic 
heating. 

3.3. Contribution to the temperature 

The question of whether the viscous heating has a significant 
impact on the overall temperature stratification remains. To ad¬ 
dress this question, one must also consider the plasma density at 
the places where the heating occurs in order to understand how 
much this volumetric heating can increase the local temperature. 
We show this in Fig.[^by looking at the Ohmic and viscous heat¬ 
ing rates per particle. We And that the mean Ohmic heating per 
particle (red) dominates over the viscous heating (blue) already 
for the full FOV (dashed lines) from a height of 8 Mm upwards. 

When restricting the FOV to the hot AR core loops SL 1- 
3 (solid lines in Fig.[^, we see a significant increase of the 

Article number, page 3 of[^ 
































A&A proofs: manuscript no. Paper 


Ohmic heating per particle in the heights between 5 and 40 Mm, 
where these loops are located. We also see a slight increase of 
the viscous heating per particle at these heights, but in the up¬ 
per corona the viscous heating per particle remains roughly un¬ 
changed. Above the smaller hot AR core FOV, the Ohmic heat¬ 
ing clearly dominates over the viscous heating by a factor of 
10 and more up to a height of roughly 30 Mm. Furthermore, in 
this region the mean of the Ohmic heating per particle (red solid 
line) also reaches over any local maximum of the viscous heating 
(blue dotted line). In the upper corona above 50 Mm, we see the 
same tendency as in Fig.[^ namely that the Ohmic heating per 
particle significantly drops below the viscous heating for the re¬ 
stricted FOVs, while this is not the case for the full FOV (dashed 
lines). 

The viscous heating can gain dominance over the Ohmic 
heating above hot closed-field regions that are strongly heated 
by Ohmic dissipation. We see that the viscous heating is slightly 
increased above the crossing point of Ohmic and viscous heat¬ 
ing compared to the full box FOV, while it is slightly decreased 
below this crossing point. In contrast to the restricted FOV, the 
global mean viscous heating profile has a rather constant scale 
height above 20 Mm. This behavior of the viscous heating could 
be explained by downflows from the upper corona that are dis¬ 
turbed by the magnetically closed AR core. Above 80 Mm the 
height profiles of the Ohmic heating per particle dominate again 
over those of the viscous heating with a factor larger than 2. 

Nonetheless, the maximum temperature stratification re¬ 
mains mainly constant in the corona above 20 Mm. The rapid 
rise of this temperature profile coincides exactly with the strong 
rise in Ohmic and viscous heating in the lower 3-7 Mm. The 
strong gradient of the temperature ends exactly where the strong 
increase of the viscous heating ends and does not follow the ad¬ 
ditional rise in the Ohmic heating towards 10 Mm. 

This overview indicates that the main contribution to the high 
coronal plasma temperatures is by the Ohmic heating, while the 
viscous heating seems to contribute mainly in higher density ar¬ 
eas where the local input to the internal energy would not result 
in such a high temperature rise. Still, the location and shape of 
the transition region temperature rise seems to be determined by 
the Ohmic and viscous heating equally, where both of them de¬ 
liver similar amounts of internal energy. 

We And a peak in viscous heating per particle at the base of 
the transition region (around 6 Mm), where coronal downflows 
(of approximately 200 km/s) are practically stopped at the top of 
the chromosphere by momentum conservation. The strong ve¬ 
locity gradients on the grid scale are numerically stable over 
the whole simulation run. These downflows are caused partly 
by cooling material converting potential into kinetic energy. In 
addition, some material that was previously accelerated by the 
Ohmic heat input leads to an increase in the plasma pressure and 
hence pushes material via adiabatic expansion. Therefore, some 
of the later viscous heating is ultimately caused by the initial 
Ohmic dissipation. 

The location of the strongest Ohmic heating per particle 
along individual loops can be inferred from Fig.[^ where the red 
solid line represents the average over the volume containing the 
three smaller loops SL 1-3 that show the strongest heating rates 
at a height of 8 Mm and have apex heights of 11-13 Mm. Fur¬ 
thermore, the hot and larger loop CL 1 has its strongest heating 
per particle between 6 and 7 Mm and reaches an apex height of 
17 Mm, while the much cooler and less EUV-emissive loop CL 2 
is 21 Mm high but its heating also peaks around 6 to 7 Mm in the 
lower corona. The heating per particle is obviously strongest at 
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some intermediate height, which we would like to refer to as the 
“legs” of the loops. 


3.4. Energy supplied into the model domain 

Perturbations in the magnetic held can travel into the corona 
with the Alfven velocity. Depending on the speed of the driv¬ 
ing motions, these perturbations will follow the magnetic held 
as transversal or longitudinal waves that could be dissipated in 
the solar atmosphere, which is referred to as “AC” heating. If 
the driver only gradually changes a held line, currents can be in¬ 
duced that compensate these slow perturbations, which can be 
seen as a diffusional process that belongs to “DC” heating (e.g., 
|Klimchuk|2006| . 

The field-line braiding introduced by our lower boundary 
creates the Poynting flux that powers the coronal energy input. In 
the photosphere, the applied horizontal motions shuffle the foot- 
points of magnetic held lines with velocities of about 1 km/s. 
This bends the magnetic held that starts to deviate from the ini¬ 
tially force-free state and these perturbations can then propagate 
with roughly the Alfven speed along the held lines. As a result, a 
net upwards directed Poynting flux carries magnetic energy from 
the photosphere into the corona, where it can then be dissipated. 
In the right panel of Fig.[T]we display the vertical Poynting flux 
at the base at the corona, which is located at a height of about 
3 Mm. 

The vertical component of the Poynting flux is highly struc¬ 
tured at the photospheric level, as are the magnetic held and the 
horizontal velocities. A change of the photospheric advection di¬ 
rection of the magnetic held can already cause a sign reversal in 
the vertical Poynting flux. This is why we see strong fluctuations 
(blue and red stripes in Fig.[2c) above the main polarities in the 
AR core even at the base of the corona. 

With increasing height, the magnetic held is smoothed out 
by the magnetic expansion and the contrast in the Poynting flux 
also smoothes quickly. In order to become independent of the 
large fluctuations in the photosphere, and also to have a reliable 
and smooth estimate of the magnetic held lines, we choose the 
height where the temperature reaches a value of 75000 K as point 
of reference for the footpoints of a held line. This temperature 
lies above the chromosphere but still within the transition region 
that is typically located at heights between 3 Mm and 5 Mm in 
our model. 

An observed coronal loop is typically a roughly semi¬ 
circular structure with a held inclination of 0 to 30 degrees in 
the lower 3 Mm of the atmosphere, i.e., below the transition re¬ 
gion. Because the magnetic held can vary strongly there, a good 
estimate of the loop length L would disregard a held line that ex¬ 
pands to remote regions by following directional changes in the 
chromospheric magnetic field. For this analysis, the loop length 
L is defined as the distance between both loop footpoints plus 
the distance to the photosphere, which we compute by assuming 
an average inclination of 20 degrees for the magnetic held below 
the transition region. Only then can we compare our model data 
to scaling laws that are usually based on observations and only 
use the visible coronal part of a loop to determine its length. If 
we use another average inclination between 0 and 30 degrees be¬ 
low the transition region, the results of this analysis deviate only 
marginally because the loop length is dominated by the coronal 
part. 




Ph.-A. Bourdin, S. Bingert, H. Peter: Coronal energy input and dissipation in an AR model 



Fig. 4. Vertical net Poynting flux into a field line versus the volumetric 
Ohmic heating integrated along that field line. Red data points show 
the absolute value of the Poynting flux wherever the Poynting flux is 
negative, dark blue dots show the data points where the Poynting flux 
is positive. The gray dashed line indicates the equality, where all of the 
incoming Poynting flux would be converted to Ohmic heat. We apply 
a least-absolute-deviation (LAD) fit to the positive net flux data points 
(black line). The green stars indicate the mean Ohmic heating within 
the corresponding bins in Poynting flux. See Sect. |4.1.T] 

4. Statistical analyses of coronal field lines 

For a statistical analysis of coronal field lines we first need to 
select a sample of field lines that represent the AR corona. From 
our model we extract not only the field lines along hot coronal 
loops, but also field lines rooted all over the AR region core. 
The seed points for the field-line tracing are distributed along a 
vertical cut through the domain, located roughly in the middle 
between the two main polarities of the AR. Selected field lines 
are a) longer than 18 Mm in total, b) shorter than 150 Mm, and 
c) reach at least to the given minimum temperature of 75000 K 
somewhere below 18 Mm height. If a field line does not fulfill 
one of these conditions it either closes below the corona, it is 
located in a cool quiet region and would not produce coronal 
emission in the spectral lines we are interested in, or it spans 
high up into the corona where the magnetic perturbations from 
the bottom layer might not yet have sufficiently propagated. We 
also exclude open field lines from the analysis. In addition to this 
rather “global” ensemble, for each hot AR core loop we added 
one field line that crosses the maximum intensity of that loop, as 
well as about 200 field lines neighboring these loops (see Fig.[^. 
The 67000 field lines that we use for this analysis have apex 
heights of about 6 to 80 Mm, while the most prominent loops 
seen in synthesized EUV emission reach heights between 10 and 
22 Mm. 

4. 1. Coronal energy conversion 

In Sect.[^ we analyzed the coronal energy conversion process 
only with spatial averages over different subvolumes. In this sec¬ 
tion we now turn to the coronal energy conversion process, in 



Fig. 5. Same as Fig.|^ except that we show here the mean Ohmic heat¬ 
ing per particle along a field line. 


particular the Ohmic heating and its relation to the photospheric 
driving and to coronal temperatures, by looking at ensembles of 
coronal field lines. 

4.1.1. Poynting flux as energy source 

We investigate the influence of the Poynting flux in the transi¬ 
tion region on the coronal energy input. The Poynting flux into a 
field line is computed as the sum of the directed flux at both field¬ 
line footpoints, while an upward flux is positive and a downward 
flux is negative. This sum gives a measure for the net vertical 
Poynting flux that actually remains in the corona, if that value 
is positive (58000 field lines). For the one snapshot we use for 
this study, about one seventh of the selected field lines show a 
negative net flux (9000 field lines), which stands for a magnetic 
energy outflow from the corona downwards. The sign of the ver¬ 
tical Poynting flux will fluctuate over time, depending on the 
evolution of the photospheric driving motions. Nonetheless, also 
for a single snapshot, the horizontal average of the divergence of 
the Poynting flux is positive at all heights within the hot AR core 
loops FOV, see Fig.[^ 

In Fig.|^we relate the vertical net Poynting flux into a field 
line to the volumetric Ohmic heating integrated between both 
field-line footpoints, not taking into account the part below the 
transition region. The main bulk of the population is of positive 
net Poynting flux (dark blue). Negative net Poynting flux also 
occurs (red), but the strongest Poynting fluxes are all positive, as 
can be seen at the right end of the distribution in Fig.[^ The large 
scatter is due to the evolution of the photospheric driving where 
we look only at one snapshot in time. A negative net Poynting 
flux (red points in Fig.[^ can be explained, for example, by coro¬ 
nal induction or Lorentz work that increases the magnetic en¬ 
ergy locally, which can then be transported downwards. Another 
likely process is that the negative flux (produced at the lower 
boundary) cancels out with positive flux from neighboring field 
lines due to the magnetic expansion in the corona. In total, the 
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mean Poynting flux for the full FOV of our model is nonetheless 
positive at every height. 

All hot AR core loops (marked by symbols) lie within the 
positive net flux population on the upper right except for SL3, 
which has a negative Poynting flux, see Sect. 4.1.3| 

We fit a scaling law with an exponent of a = 0.47+0.06 to the 
positive data points (see Fig.|^. The error is computed from the 
standard deviation of the LAD fit and by estimating the width of 
the data distribution to four orders of magnitude in the horizontal 
direction. While the scatter of the individual data points is quite 
large, the trend found by the fit is still clear. This is supported by 
grouping the data points into bins of the vertical Poynting flux 
(x-axis in Figs.andand calculating the mean of the Ohmic 
heating in these bins. These data points are shown by green stars 
in Figs. 0 and 1^ This shows a small scatter and a clear trend 
consistent with the power-law fit. The positive slope of the fit is 
consistent with a relation between the volumetric Ohmic heating 
integrated along a coronal field line and the vertical Poynting 
flux at the field-line footpoints of the form 


Fohm ~ 


All in all, this is good because positive exponents a and fi support 
a possible correlation between the net Poynting flux into a field 
line and its temperature, via the Ohmic heating in the corona, 
which we would expect. 

We get a similar result when we use the Ohmic heating per 
particle instead of the volumetric heating for this analysis. 

4.1.2. Ohmic heating integrated along field lines 


For the same ensemble of field lines as in Sect. |4. l.T] we inves¬ 
tigate in Fig.|^the Ohmic heating per particle. We find here that 
the population of hot field lines mainly follows a scaling law of 
the form 


T ^ 

^ max ^ Ohm ? 


( 2 ) 


where the maximum temperature T^ax relates to the Ohmic heat¬ 
ing integrated along that field line Fohm that we calculate from 
the volumetric Ohmic heating Hohm as: 


( 1 ) 


Fohm ~ 


where !P = 77(VxB)xB + (uxB)xB//i with the velocity u, the 
magnetic field vector B, the magnetic diffusivity rj, and the per¬ 
meability fi. Above a height of 2 Mm, where all the footpoints of 
our field lines are located, the Poynting flux is largely dominated 
by the second term that reflects transversal perturbations in the 
magnetic field. As the main distribution of data points is clearly 
located below the equality line, we conclude that not all of the 
Poynting flux in our model is usually converted to Ohmic heat. 
The hot AR core loops SL 1-3 and CL 2 convert a larger part of 
their incoming Poynting flux into Ohmic heat than the average 
field lines do. 

Generally, the Ohmic heating integrated along field lines is 
smaller than the Poynting flux, although the energy is conserved 
because not all Poynting flux needs to be converted into Ohmic 
heat (or dissipated at all). Some of the remaining upwards mag¬ 
netic energy transport that is not dissipated can leave the sim¬ 
ulation box though the upper boundary or it can be converted 
into magnetic field stress that can let the coronal magnetic field 
evolve away from a potential field configuration, or it can lead to 
a rise of magnetic structures into the upper corona. For instance, 
we found that the large AR core loops CL 1 and 2 are still rising 
with about 2km/s at the loop top ( |Bourdin et al.|2013] l. Accord¬ 
ingly, we find a net upwards mass transport in the middle of the 
AR core between the two main polarities. 

We now perform the same analysis with the mean Ohmic 
heating per particle instead of the volumetric heating, see Fig.[^ 
Again, we overplot a LAD fit to a scaling law that gives an ex¬ 
ponent of y6 = 0.38 + 0.07 (black dotted line), which indicates 
a weaker correlation. We also find a 14% higher variance in the 
logarithmic data, and hence we see more scatter in Fig.j^than in 
Fig.|^because the heating per particle also depends on the par¬ 
ticle density n, which plays a role when relating the magnetic 
energy input not only to the dissipated magnetic energy, but also 
to the actual heating of the coronal plasma. 

For the strongest heated loops SL 1 +2 we find a heating per 
particle of about 10“^ eV /s, while the average for field lines with 
the same magnetic energy input is about 10 times smaller (see 
the green stars around 10“^eV/s below the black symbols in 
Fig.|}. 

n both Figs. 1^ an d0 the hot loops (symbols) clearly have a 
stronger heating than the scaling law fitted through all data points 
would predict (black line) and also a higher temperature than the 
average of field lines with a similar Poynting influx (green stars). 
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The exponent of 2/7 is predicted by the RTV scaling law (named 
after the authors of |Rosner etal.|1978] l. 

While the distribution of data points in Fig.[^ is relatively 
narrow for the high heating rates and high temperatures, the dis¬ 
tribution becomes much broader for field lines with a maximum 
temperature below 1 MK. These cooler field lines have insuffi¬ 
cient heating to become hot and are not quasi-static. Therefore, 
these will not follow the RTV scaling law and one expects scat¬ 
ter. 

However, the hot end of the population of field lines with 
more than 100 W/m^ of integrated Ohmic heating follows the 
slope of 2/7 and is comparable with a scaling law, see upper left 
bound in the right half of Fig.which represents the hot coronal 
loops that the RTV scaling was originally designed for. 

Still, the large scatter in Fig.|^shows that the scaling laws of 
|Rosner etal.| ( |1978) l can only be considered an average variation. 
Any realistic model considering the temporal variability and the 
large range of heating rates will produce loops that do not fol¬ 
low the RTV laws. In this sense the RTV scaling laws should be 
considered only on a statistical basis. 

4.1.3. Mean Ohmic heating 

Now we investigate the relation between the average heating per 
particle and temperature along a field line. For that, we need to 
replace the volumetric Ohmic heating Ffohm in Eqn.[^ with the 
Ohmic heating per particle Hohmin- The RTV scaling law to¬ 
gether with the ideal gas law p ~ nkgT predicts a slope of 2/3: 


{T) ~ {Hohmlnf^ 


(4) 


The relation in Fig. [7] between the mean Ohmic heating per 
particle and the mean temperature along a field line is clearly 
visible for field-line lengths between 60 and 110 Mm (turquoise 
to red colors), which follow nicely the predicted slope of 2/3 
(black dotted line). For field-line lengths below 55 Mm, as well 
as for mean temperatures below 0.2 MK, we find large scatter. 

Nonetheless, here we see a clearer match to the predicted 
slope, even for groups of different field-line lengths above 
60 Mm (see red, yellow, green and turquoise populations in 
Fig.|7]), than in Fig.|^ This shows again, as discussed in Fig.|^ 
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Fig. 6. Volumetric Ohmic heating integrated along a field line (Fohm) 
versus its maximum Temperature The hottest loops (marked 

with symbols) all have a high Ohmic heating. The three loops with the 
strongest EUV emission (SL 1+2 and CL 1) lie exactly at the local max¬ 
ima of their surrounding population. The dashed black line indicates the 
exponent of 2/7 with the parameters of the loop SL 1 as anchor point, 
see Sect. l4.L2l 


Fig. 7. Mean Ohmic heating per particle along a field line {{Hohmin)) 
versus its mean Temperature {{T)). The mean is computed here as the 
heating integrated along a field line divided by the integration path 
length. The field lines with a heating rate above 10“* eV/s are roughly 
ordered according to their length. The black dotted line indicates a slope 
of 2/3, see Sect. |4.L^ 


that the density is a relevant ingredient in order to describe not 
only the coronal heating, but also the resulting plasma tempera¬ 
ture along a field line. 

The short and hot loops SL 1 +2 show the strongest mean 
heating, while the longer loops CL 1+2 have a 10 times lower 
average heating, even though their lengths (35 and 70 Mm) dif¬ 
fer only by a factor of 2. Still, they have roughly the same aver¬ 
age temperature, because the very efficient heat conduction and 
the radiative losses act like a thermostat. We see the same results 
when we look at the volumetric heating (Fig.|^ instead of the 
heating per particle. 

The loop SL3 is relatively cool and does not show much 
emission in EUV and X-rays and is therefore not a typical hot 
coronal AR loop (cf. [Bourdin et al.||2OT3] l. The lower mean 
Ohmic heating can be explained because SL 3 is the only loop 
with a negative net Poynting flux, see Sect. [4.1.1 CL 2 is not a 
developed hot loop and lies outside the main coronal emission 
structures. 


CL 1 lies farther away from the equality line (gray line in 
Fig.[^, which indicates that this loop uses a smaller fraction of 
the incoming Poynting flux for Ohmic heating than the other hot 
AR core loops. Owing to the particularly high incoming Poynt¬ 
ing flux, CL 1 still can reach Ohmic heating rates as high as the 
other loops. 


4.1.4. Relation to magnetic field strength 

Also of high interest is the connection between the field-line 
length and the total Ohmic heat input (which sets a field line’s 
tendency towards high or low coronal temperature) in relation 
to the magnetic flux density at the footpoints of the field line. 
As mentioned in Sect.[2.1[ there are strong fluctuations of the 


vertical magnetic flux density in the photospheric layer. It is im¬ 
portant that this analysis is independent of these strong fluctua¬ 
tions that would anyway not reach into the corona. Therefore we 
choose 300 km as the reference height for the magnetic flux into 
a field line, which is roughly one photospheric pressure scale- 
height above the photosphere, but still below the chromosphere 
in our model. 

In Fig.[^we display the mean of the vertical magnetic flux 
density values at both field-line footpoints versus the volumet¬ 
ric Ohmic heating integrated along each field line. The field¬ 
line length (color coded) is taken between the reference heights 
where the temperature reaches 75000 K, plus the extension to¬ 
wards the photosphere. 

Below a critical magnetic flux density of roughly 200 G (gray 
dashed line in Fig.[^, we find a wide distribution of flux densities 
with no general correlation to the integrated Ohmic heating or 
to the field-line length. Nevertheless, almost all strongly heated 
field lines are located above this critical flux density. We do not 
see a strong correlation between the footpoint flux density and 
the coronal heating rate, but the field-line lengths clearly separate 
and range from the long field lines (100 Mm, red at the lower 
right) with relatively low Ohmic heating to the short field lines 
(30 Mm, dark blue at the upper left) with the strongest Ohmic 
heating. 

Taken all together, this supports the interpretation that the 
magnetic energy input into a field line is responsible for the coro¬ 
nal Ohmic heating and that a threshold of about 200 G in the 
footpoint flux density is a prerequisite to providing a sufficient 
net vertical Poynting flux. A sufficiently strong flux density also 
allows for a field line that reaches into the corona, instead of 
closing at low height to one of the many neighboring opposite 
polarities of less than 200 G. 
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Fig. 8. Magnetic flux density at the footpoints located at z=300km 
height versus volumetric Ohmic heating integrated along the field line. 
The hot loops SL 1+2 lie within the high heating rates population. The 
vertical gray dashed line indicates the magnetic flux density threshold 
of 200 G, see Sect.|4T4] 


5. Conclusions 

We found the low density to be a prerequisite for a strong rise 
of the Ohmic heating per particle and hence the coronal loop-top 
temperature (see Fig.[^. This is interesting because a lower pres¬ 
sure could instead be the result of strong heating that expands the 
plasma adiabatically. However, we see equally strong volumet¬ 
ric Ohmic heating in large parts of our model corona and only 
where the density is low does the temperature rise later. This 
supports the idea that in the real corona the volumetric heating 
is also distributed in space more evenly than the observations of 
very confined regions with EUV-bright loops would suggest. 

The RTV temperature scaling law fits roughly to our hottest 
and strongest heated coronal field lines, while this scaling law 
does not hold for less heated or cooler field lines. Looking only 
at the field lines originating at the center of the coronal loops, 
our model data deviates from the RTV temperature scaling law 
by 25 %, see symbols in Fig.|^ This indicates that the RTV scal¬ 
ing law is not per se applicable to any coronal field line, but 
is roughly consistent with EUV-emissive coronal loops that this 
scaling law was originally fitted to. 

Regarding the discussion whether loops are heated from their 
footpoints or along the loop, we can say that in our model loops, 
the volumetric heating is of course stratified following the ex¬ 
pansion of the magnetic field (see Fig.|^. However, when look¬ 
ing at the heating rates per particle (Fig.|^, we see that it is 
actually strongest around 9 Mm height. Therefore, speaking of 
“footpoint heating” (in absolute volumetric values) is physically 
correct, while at the same time it is irrelevant for basic properties 
of the coronal loops, instead it would be better to better speak 
of “loop-leg heating”, because even the term “loop-top heating” 
would be incorrect here (cf. [Priest etal. 11998) 1. 

This finding is perfectly consistent with the field-line braid¬ 
ing mechanism, because the photospheric horizontal motions 
will have the strongest impact through currents generated by 
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deformation or displacement of a magnetic field line where the 
field is still rather vertical. The resulting heating per particle is 
strongest where the density is low enough. Both conditions are 
optimally fulfilled for coronal loops in their “legs”, which we 
refer to as the lower coronal part of a loop at heights around 
5-lOMm. 

Fundamental in this model is that the heating is due to a slow 
magnetic energy dissipation process, which would be described 
best as steady Ohmic heating due to quasi-static magnetic field 
changes induced by photospheric advection leading to a slow 
reconfiguration with simultaneously generated electric currents. 
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